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Abstract: We construct explicitly time-dependent exact solutions to the field equations of 
6D gauged chiral supergravity, compactified to 4D in the presence of up to two 3-branes 
situated within the extra dimensions. The solutions we find are scaling solutions, and are 
plausibly attractors which represent the late-time evolution of a broad class of initial con- 
ditions. By matching their near-brane boundary conditions to physical brane properties we 
argue that these solutions (together with the known maximally-symmetric solutions and a 
new class of non-Lorentz-invariant static solutions, which we also present here) describe the 
bulk geometry between a pair of 3-branes with non-trivial on-brane equations of state. 
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1. Introduction 

The discovery of the existence of D-branes [Q within string theory has led to a fundamental 
rethinking of the kinds of effective theories which can describe the low-energy limit of a 
fundamental theory. The broadening of the collective mind to which this rethinking has led 
has permitted new progress to be made on a number of the 'naturalness' issues which seem to 
plague our understanding of the low-energy theories — i.e. the Standard Model plus General 
Relativity — which describe the world we see around us. 

This progress has come about largely because the possibility of trapping low-energy 
particles on surfaces within higher-dimensional spacetime changes how one must think about 
naturalness issues. For instance, an important property of 4 dimensions is the equivalence of 
a vacuum energy with a cosmological constant, and so also with nonzero 4D curvature. This 
connection underlies the cosmological constant problem |5| , which amounts to the difficulty 
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in understanding why the observed universe can be so flat given that quantum fluctuations 
generically make the vacuum energy enormously large. 

This connection can be broken in higher-dimensional brane configurations, inasmuch 
as known higher-dimensional solutions show that large 4D energy sources can coexist with 
zero 4D curvature, as has been pointed out for pure gravity with co-dimension one U 
and co-dimension two || [7], || geometries, as well as for co-dimension two geometries within 
supergravity |], 10, ll], O]. They can coexist in this way because within the extra-dimensional 



context the 4D energy density turns out to source the curvature of the extra dimensions rather 
than the curvature of the 4 dimensions within which the energy density exists. 

Whether or not this observation can lead to a solution of the cosmological constant 
problem depends on identifying the assumptions which are required in order to ensure that 
the 4D curvature is sufficiently small, and on whether or not these choices are natural or if they 
must be fine tuned. The hope is that the reformulation of the cosmological constant problem 
in this way can usefully recast the issues into a form which might allow progress that was 
not possible in the 4D context. 6D supergravity provides a particularly attractive framework 
within which to test these tuning issues in detail, since in six dimensions it is possible to have 
internal geometries which are large enough to allow the Casimir energy in these dimensions 
to be of order the size of the observed cosmological constant, allowing a variety of attractive 
phenomenological possibilities J13| ] . In this paper we use 6D chiral gauged supergravity to 
explore one of the fine-tuning issues which arise within this extra-dimensional context. 

There are two different kinds of fine-tunings against which one must be vigilant within 
these extra-dimensional models. The first of these is the requirement to tune against quantum 
fluctuations. That is, if parameters are chosen to ensure that 4 dimensions are flat, do these 
dimensions remain flat once the scale of the theory is lowered by integrating out particles 
having successively lower masses. This is the traditional cosmological constant problem, 
and this paper has nothing to add to the ongoing investigations as to whether the higher- 
dimensional models require tuning in this way |l4]| . Our focus here is instead on the second 
fine-tuning issue, which arises in any formulation (such as the higher-dimensional scenarios 
of interest) for which the Dark Energy density is time-dependent. This second question asks 
whether having an acceptable present-day cosmology requires an inordinately finely tuned 
adjustment of initial conditions before the advent of the present epoch's Big Bang cosmology. 

In order to do this we shall construct several nontrivial exact solutions of the Nishino- 
Sezgin 6D chiral gauged supergravity [15, 16, 17]. These solutions will correspond, in general, 
to the introduction of non-trivial matter sources on the brane, but in specific cases they 
describe the response of the bulk geometry to pure tension branes whose tensions are not 
otherwise fine-tuned. In this sense we are probing the cosmological constant problem in 
these six dimensional compactifications, by inferring the brane geometries induced, via the 
bulk geometry, from arbitrary brane tensions. From a cosmological point of view it is also 
interesting to have explicit solutions for branes with arbitrary equations of state, and couplings 



1 See also ref. [hq] for a class of solutions which are very similar to some of the ones we present here. 
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to the bulk fields such as the dilaton. We take a pragmatic approach of first looking for bulk 
solutions, and then inferring the properties of the branes for which these would be consistent 
solutions. 

The remainder of the paper is organized as follows: we close this section with a brief 
summary of the relevant field equations. §2 then describes the properties of a 3-parameter 
class of static solutions having two compact dimensions. They include the 2-parameter fam- 
ily of previously-known solutions which are maximally symmetric in the four noncompact 
dimensions, as well as a new set of static solutions which break the 4D Lorentz symmetry. 
These solutions generically contain singularities which we can associate with the positions of 
two source branes. §3 derives a broad class of time-dependent solutions describing geometries 
for which the sizes of the various dimensions scale with a power of time. A related class of 
solutions describing one brane within two noncompact dimensions is then described in §4. 
In §5 we argue that these solutions are rich enough to describe the late-time behaviour of a 
large class of source branes. We provide the matching conditions which allow the asymptotic 
behaviour of the bulk fields to be related to physical properties of these branes like energy 
density, pressure and dilaton coupling. We close in §6 with a brief summary and conclusions. 



Field equations 



The bosonic part of the Lagrangian density for 6D chiral gauged supergravity []15| |16|, |17|] is 
given by 2 



-g 2k 2 



Rmn + d M 4>d N (, 



K r MN r 



2g 2 



where (ft is the 6D scalar dilaton and F = dA is the appropriate field strength for the gauge 
potential, A M , which gauges a specific abelian .R-symmetry, whose gauge coupling, g, has 
dimensions of inverse mass. We keep the 6D Planck scale, k 2 = M^~ 4 , explicit for ease of 
comparison with the various conventions which are used in the literature. These expressions 
set to zero some of the bosonic fields of 6D supergravity, including possible matter hyperscalars 
and gauge potentials, <5 a = = 0, and Kalb-Ramond fields, G MNP = — a choice which 
is consistent with the corresponding field equations pi]]. 
This action leads to the following field equations: 



D M NT* F MN ) = (Maxwell) 



K 

T 



□ </>+— e" 9 ' F MN F MN - = (Dilaton) 



K 2 



Rmn + d M (f> d N <p + K 2 e"^ F MP F N P + -(□<£) g MN = (Einstein) 



(1.2) 



2 The curvature conventions used here are those of Weinberg's book [[L9[ , and differ from those of MTW 
only by an overall sign in the Riemann tensor. 
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to which the bulk of the remainder of this paper is dedicated to solving, in order to find explicit 
static and time-dependent compactifications to four dimensions. To this end we divide the 
six coordinates x M , M = 0, . . . , 5, into 4D coordinates x M , \i = 0, . . . , 3, and 2D coordinates 
x m , m = 4, 5. When required, the three spatial coordinates of the noncompact 4 dimensions 
are denoted x % , i = 1,2, 3. 

2. Static solutions 

We start by describing a broad class of static compactifications. For some of these solutions the 
noncompact 4 dimensions are maximally symmetric, and these solutions have been described 
in the literature. We supplement these with new static solutions which break the 4D Lorentz 
symmetry. 

2.1 Ansatze 

We start by writing out the field equations for configurations which are (i) time- independent; 
(ii) translation and rotation invariant in the 3 noncompact spatial dimensions; and (Hi) are 
axially symmetric in the extra dimensions. That is, we take 

ds 2 = _ e M>(„) dt 2 + e 2a( V ) § .. dx i dxj + e 2v(v) dr] 2 + ^b( V ) dQ 2 

A e = a (v) and = e^ v) , (2.1) 

leading to the following system of coupled ordinary differential equations: 

a'0 + (w' + 3 a' - v' - V - a! e = (Maxwell) 



?" +(w' + 3a'-v' + b') 0' + ^- e- 2l ^{a' e f - \ e 2 ^ = (Dilaton) 

2 2 

w" + (w' + 3 a' - v' + V)w' - y e" 2i -^{a' e f + ^ e 2i+0 = {tt Einstein) 

2 2 

a" + (w 1 + 3 of - v + V)a! - ^ e" 2 ^^ \a' e ) 2 + ^ e 2i+0 = (ii Einstein) 

b" + (w 1 + 3 a! - v' + V) b f + ^f e- 2l -v(a' e ) 2 + ^ e 2 ^ = (66 Einstein) 



w" + 3 a" + b" + (w'Y + 3 (a 1 ) 2 + (b 1 ) 2 + (0'f ( W Einstein) 

-{w' + 3 a' + V) v' + ^f e- 2b -£(a' e ) 2 + ^ e 2 ^ = . (2.2) 

Here ' denotes a derivative with respect to r\. These equations must be supplemented with 
boundary conditions at the locations of the two branes. 

Although this appears to provide 6 equations for the 6 unknown functions w, v, b, (p 
and a$, this is deceptive because we can ensure that one of these functions (say v) takes any 
particular form simply by appropriately changing the coordinate r]. However, one combination 
of these 5 equations, found by taking the combination (7777) — (tt) — 3(ii) — (66), can be 
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thought of as a 'constraint' on the evolution of the fields into the r\ direction because all 
second derivatives, (d/drj) 2 , drop out 3 : 



k 2 e- 2b -Ha' ft ? 



6 {a' f + 6 a'w' + 6 a'b' + 2 b'w' - (fif + 



f\2 



(2.3) 



As a consequence of the Bianchi identities, one of the remaining 5 equations is then not inde- 
pendent and can be derived from the derivative of the constraint and the other 4 equations. 

A considerable simplification to this system of equations can be obtained by performing 
the redefinitions 

w = 3£ + (y - x)/4 

a = -€+(y-x)/4 
v = InN + (5y - x + 2z) /4 
b = (3x + y + 2z)/4 
= (x-y-2z)/2. 



Then one can demonstrate that the system of equations ( |2.2[ ) and the constraint (^ 
from the Euler-Lagrange equations associated with the action 



S 



ch] 



N- 1 [(x'f - (y') 2 + (z'f + 12(£') 2 + K 2 e- 2x (a' e ) 2 ] + 



N 



V 

K 2 



(2.4) 
follow 

(2.5) 



here plays the role of a Lagrange multiplier and can be set to unity after variation. This 
system can be solved exactly, giving (with N = 1): 



061 



q J dr/e 
nq 



2.t 



cosh[Ai (rj - rji)] 



cosh[A 2 [r] - r] 2 )} 



k\ 2 
z = z + A 3 rj 

£ = Co + ^4 T) 



(2.6) 



and where the constraint equation amounts to the condition A2 = Af + X 2 + 12 X 2 . We can set 
two of the parameters 771 , 772, Zq, £0 to zero by coordinate rescalings without loss of generality. 

2.2 Lorentz-invariant solutions 

A special case of these solutions is A4 = 0, for which £' = and so a' = w' . The resulting 
solutions are therefore Lorentz invariant, and reduce to those that have been found earlier in 



3 This equation follows directly from writing the (7777) Einstein equation in terms of the Einstein tensor, 



0, and it is the Bianchi identity which ensures that it holds for all rj once it is imposed on 



'initial conditions' at 77 = 770- This is the Hamiltonian constraint associated with the lapse function expTj in 
the spacelike ADM formalism of gravity. 



- 5 - 



the literature g|, [To], H H- 4 In thi s 

case the general solutions can be written as 
ds 2 = W 2 (i]) Tj/ju/ dx^dx u + A 2 (rj) [w 8 (r?) drf + d# 2 

F v e = (^) e- x ^ and e"* = W 2 e A ^ , (2.7) 



where 



W 4 



/ K 2 q\ 2 \ cosh[Ai(r/ - 771)] 
V 2g\ x J cosh[A 2 (r? - 772)] 



-4~ 4 = (^l^^e-^cosh^Ai^-^^cosh^lT?-^], (2.8) 

along with the constraint A 2 . = A 2 + A3. 

A further special case is obtained for A3 = A4 = 0, and r\\ = 772 for which <p = ipo and 
w = a = wq are constants. Changing variables to proper distance, dp = e^'drj, the form 
field becomes 

F pe = a' g = ±4^°B(p), (2.9) 
where primes now indicate differentiation with respect to p and B(p) = e b ^ satisfies 



^l = b" + (b') 2 = -^-e^, (2.10) 



with solution B(p) = Bosm^ge^^p/nj. These are the rugby-ball generalizations ||] of the 
older Salam-Sezgin solution [15] to gauged chiral 6D supergravity. 

2.3 Asymptotic Forms 

These solutions describe geometries which become singular at r/ = ±00, which are interpreted 
as being the positions of the 3-branes which source this configuration. Since it is ultimately the 
internal structure of the brane (if any) which is responsible for resolving these singularities, 
one expects this structure to be related to the asymptotic limit of the above solutions in 
the near-brane limit. We therefore pause here to outline what this asymptotic near-brane 
behaviour is. 

To this end it is useful to adopt Gaussian-normal (GN) coordinates near the brane for 
which ds 2 = g a b dx a dx b + dp 2 , where x a denotes the 5 other coordinates, {x a } = {x M ,6>} = 
{t,x l ,6}, with the brane position being described by p = 0. We take the asymptotic form of 
the bulk fields in the near-brane limit (p ^ £) to be generically given by a power law [23 

ds 2 ~ -[c w (H lP r] 2 dt 2 + [c a {H lP ) a ] 2 Sij dx'dx' + dp 2 + [c e {H 1 p)^ l ] 2 p 2 de 2 

~ c^Hxpf and F pe ~ c f (H lP y , (2.11) 



4 The conventions of ref. [£2| may be obtained from ours by taking Rain — > —Rmn, 4> —> ~ 0/2 and n 2 = 1 /2, 
while those of pd| differ from those here only by the choice k 2 = 1. 
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where to, a, f3, p, 7, c w , c a , ce, and c/ are constants, and H\ is an arbitrary scale. With 
these choices the extrinsic curvature, = ^ d p g a b, of the constant-/? surfaces becomes 

iCt**-, K^^-S'j and (2.12) 
P P P 

up to contributions that are subleading for small p. 

Only two of the five powers a, f3, uj, 7 and p defined above are independent, since the 
bulk field equations impose the following two conditions amongst them 

u 2 + 3a 2 + (3 2 +p 2 = uj + 3a + (3 = 1 and 7 = p - 1 . (2.13) 

For example, explicit calculation with the general static solutions given above gives 

A 2 -Ai±4A 4 A 2 -Ai=f12A 4 , 01 .. 

5A 2 - Ai =F 2A 3 5A 2 - Ai =F 2A 3 

A 2 + 3Ai t 2A 3 , 2(A 2 -A lT 2A 3 ) 
5A 2 - Ai =F 2A 3 5A 2 - Ai =F 2A 3 

at the brane positions r] — > ±00. As is easily verified, these satisfy the expressions ( |2.13 ) 
above. 

In later sections we relate these powers to the physical properties of the branes which 
source these geometries, following arguments which generalize those of refs. [12, 0, 23]. 



3. Scaling solutions 

We now generalize the previous discussion to a new class of time-dependent scaling configu- 
rations which provide exact solutions to the same 6D field equations. The idea behind the 
construction is to assume a metric of the general form 

ds 2 = -t c N 2 dt 2 + qij (t c/2 /?Mt + i c/2+ W) (t c l 2 (3 3 dt + t c ' 2+1 dx^ , (3.1) 

for arbitrary real c, where N and (3 l represent the lapse and shift functions just as in the usual 
ADM decomposition and q^j is the metric in the n-dimensional hypersurface with normal 
vector (3 l . These quantities are assumed to be independent of time: N = N(x l ), ft = 
ft(x l ) and 

Qij — Qiji^ 1 )- This is the most general metric, up to coordinate transformations, 
preserving a single time- like homothetic Killing vector that acts as t — > Ai, gMN — ► ^ 2+c 9mn- 
The key point to recognize is that with this choice the components of the Ricci tensor 
also scale as a simple power of t: Ru oc t~ 2 , Ru oc t" 1 and Rij oc t°. Because of this, and of 
the scale-invariance of the supergravity equations, it is also possible to scale the other fields in 
the problem with t in such a way as to ensure that the field equations are also proportional to 
specific powers of t. Once this is done, then the problem of solving the field equations becomes 
an exercise in determining the profiles of the metric functions in one lower dimension. 
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The most general form for the dilaton and the gauge field necessary for eq. (3.1) to be a 
consistent solution of the equations of motion is determined by the scaling symmetry to be 



(j) = - (2 + c) log t (3.2) 

On substitution of these forms into the 6d equations of motion we find that the time depen- 
dence completely decouples as expected from the scaling symmetry, and we are left with an 
effective 5 dimensional theory for the profiles qij, fli, N, (j), Aq, A{. In practice the resulting 5d 
system is formidable and so in what follows we shall consider special cases of these scaling 
solutions. In particular we concentrate on the case where the variables only depend on one of 
the spatial dimensions (namely the radial direction 77). Furthermore we will work in the spe- 
cial case where only one component of the gauge field does not vanish, in particular A§ = 0. 
All these assumptions greatly simplify this effective theory. 



3.1 Warped scaling solutions 



From now on we shall only consider scaling solutions which are warped on a single extra 
dimension. We consider the special case where the metric is 



ds 2 = (H t) c \-e 2w ^dt 2 + e 2a ^ Sy dx'dat] + (H t) 2+C \e 2v ^dr] 2 + e 2b ^d9 2 



while 



and 



a <p(v) 



(3.4) 



(3.5) 



(H t) 2+C " 

Although this metric is not precisely of the form of eq. (|3.l| ) it can easily be shown to 
be equivalent by redefining x l = tx l . From now on we shall work with the form (|3.4j ) for 
convenience. Here Hq is a constant of dimension inverse time and the functions w, a, v, b, ag 
and ip are to be determined by solving the field equations. In some cases it will be natural to 
take Hq < so that the direction of increasing time corresponds to t — > 0~ . 

By virtue of the way the Ricci tensor scales with t, with the above ansatz all of the field 
equations reduce to the following set of coupled ordinary differential equations which govern 
the 77-dependence of the various undetermined functions. The Maxwell equation is 



a'g + w' + 3a'- b' 



tp')a' e 



0. 



while the Dilaton equation similarly becomes 



" +(w' + 3a'-v > + b') tf/ + (2 + c)(2c + 1) e 2{v ~ w) H^ 



(3.6) 



(3.7) 



2g 2 



2v+<p 
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The Ricci tensor for this class of metrics is easily computed and leads to the following com- 
ponents for the Einstein equations. The (trj) component is 

(2c + l)w' + 3a' + (2 + c)(p' = 0, (3.8) 

while the (tt) equation is 

w" + (^w'+^a'-v'+b'^w'+^e-^ia'gf + ^e 2 ^-^ 2 + | c + 4) e~ 2w+2v H 2 = . (3.9) 
The (99) equation is 

6 // +(u;'+3a / -t; , +6 , )6 , + ^e- 2fe ^(4) 2 + ^ e 2 ^-i( C +2)(2 C +l) e ^ +2 ^ 2 = 0, (3.10) 
the (ii) Einstein equation becomes 

a" + (w' + 3a' -v' + b')a' - ^e' 2 ^ (a' e ) 2 + 4 e 2v+v - \c(2c + l)e~ 2w+2v H 2 = , (3.11) 
and finally the (7777) equation is 

w" + 3a" + b" + (w'f + 3(a') 2 + (b'f + (<p') 2 - (w' + 3a' + 6') v' (3.12) 

+ 3^ e -26-^ (a / )2 + g e 2, +v _ 1 (c + 2)(2c + = Q ( 

where ' = djdr\. Equation counting proceeds much as for the static solutions, although with 
an important difference. The difference is the inclusion of time dependence, which makes 
two of the components of the Bianchi identity nontrivial rather than one. This implies that 
in this case two of the field equations are not independent of the others, rather than just 
one. Related to this is the existence in this case of two constraint equations which do not 
involve d 2 /dr/ 2 , which the Bianchi identities ensure are preserved when evolved using the field 
equations in the rj direction. These constraints can be taken to be the (trj) Einstein equation, 
eq. (|3.8j) , and the combination (7777) — (tt) — 3(H) — (99): 

(^) 2 -6(7«'+a'+6 / )a , -26'7i; , +K 2 e - 26 -^(^) 2 -^e 2l ' + ^+4( C 2 + C +l)e- 2u ' +2 ^ 2 = 0. (3.13) 

The logic to solving these equations is to use eq. (^) to solve for a and eq. ( |3.13| ) to solve 
for ag. Once this is done, we may ignore the Maxwell and (ii) Einstein equations as they are 
redundant. We are then left with three independent equations — i.e. the Dilaton and the 
(tt) and (99) Einstein equations — for evolving the remaining three independent functions 
— w, b and cp — into the 77 direction. 

As for the static case the system of equations is simplified by making the choice 

w = 3£ + (y - x)/4 

a = -f + (i/-a?)/4 

v = In N + (by - x + 2z) /4 

b = (3x + y + 2z)/A 

ip = (x-y-2z)/2. (3.14) 
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The (trj) equation now amounts to the condition 



(3.15) 



up to an irrelevant integration constant. The remaining equations of motion and (7777) — (tt) 
3(ii) — (69) follow from the Euler-Lagrange equations of the action 



S = j dr/ jiV" 1 



(x'f - (y'f + 4(1 \ C / (z'f + Je^(a' e f 



+ N 



%e 2 y-AHi(l + c + c 2 )e 



2\ p 2y-2z/c 



(3.16) 



where N plays the role of a Lagrange multiplier and can be set to unity after variation. As 
before ag and x decouple and can be explicitly integrated: 



a e = q J drj e 2x 

.-x K Q 



Ai 



cosh[Ai (77 - 771)] , 



(3.17) 



but y and z do not. Their equations follow from the reduced action 



S 



{y'f + 



+ N 



4(1 + C + C 2 ) 



/\2 



3c 2 



-\\ + \e 2 y -4H 2 (l + c + c 2 )e 



2\ p 2y-2z/c 



(3.18) 



With the choice N = 1 the resulting field equations are 



y" + %e 2 y-4H 2 (l + c + c 2 )e 



2\ p 2y-2z/e 







3cH%e 2y - 2z/c = 0. 



(3.19) 
(3.20) 



The asymptotic form of these solutions as r\ — ► ±00 is y — > A 2 r\ and z — > A 3 77 and the 
constraint implies (A^) 2 = Af + 4(1 + c + c 2 )(A^) 2 /(3c 2 ). Although we have not found closed- 
form analytic solutions of these equations, they are straightforward to integrate numerically, 
with the result having qualitatively the same form as the de Sitter solutions considered in 
[p3|1. 5 This relation with the de Sitter solutions is not surprising since these solutions may be 
obtained as the special case c = —2 of the above scaling solutions. 

3.2 Useful special cases 

There are several special cases of the previous solutions which are of particular interest. 



3 See also [M for de Sitter solutions in the non-supersymmetric case. 
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Connection to 4D scaling solutions 

If c = —1 then the metric and dilaton have the time dependence 

ds 2 = ^g fiu (y)dx»dx» + tg mn (y)dy m dy n and = ^ , (3.21) 

which implies in particular that \J—g g^ v is independent of t. 

Such a scaling solution has a simple interpretation in the limit where the extra dimensions 
are large enough to justify a description in terms of an appropriate low-energy 4D effective 
theory p§|| . In the classical limit this theory contains two massless modes, corresponding to 
the 4D metric and one combination of the dilaton, eft, and radius, r, of the extra dimensions 
(for which r 2 e^ is fixed) which parameterizes a flat direction of the 4D scalar potential. The 
above scaling solution describes a time-dependent scaling along this flat direction with a fixed 
metric in the 4D Einstein frame. (More general choices for c also rescale the Einstein-frame 
4D metric.) 

Explicit scaling solutions to the field equations of the effective 4D theory describing these 
modes are known, many of which are attractor solutions to which a broad class of initial 
conditions are drawn pq ] . The solutions found here show how to extend those of the effective 
4D theory to see the profiles of the other nonzero KK modes. Because the 4D solutions are 
attractors for the 4D field equations, we might also expect that the same may be true for the 
higher-dimensional solutions found here. 

Pure tension branes 

It is known that even for the special case of pure tension branes, maximally-symmetric solu- 
tions to the 6D field equations only exist when the tensions of the two branes are adjusted 
relative to one another |]l0| , |22| , |i~2| . One might hope that the above scaling solutions might 
describe the late time behaviour of the solutions in the case that the brane tensions are not 
adjusted in the appropriate way. 

We now show that a subset of the solutions found above can indeed describe this situation. 
In order to do so we must identify when the asymptotic form of the solutions near the branes 
have the pure tension form where a ~ w i.e. £ — > 0. Since the (ijt) implies £' = (2 + c)z'/(6c) 
there are only two circumstances for which £' — ► near a brane: 



(i) z?->0, or (n) c = -2 . (3.22) 



Special case c = oo 



Case (i) of ( 3.22Q corresponds to the special case where the geometry near the brane has a 



conical singularity, since this always requires z' — > in the near-brane limit. While this is 
always possible for one brane, in general it is not possible for both branes since the equation 
of motion for z is 

z" - 3cH 2 e 2y ~ 2z/c = 0, (3.23) 
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and so \z"\ > 0. There is however, one special case for which this can be achieved. Redefining 
= \ c \~ 1 Hq an d taking the limit c — > oo we recognize that the originally coupled system 
for y, z decouples, 

y " + 1^ _ 4 ff 2 e 2 ^ = z" = 0. (3.24) 
For which the conical solution is 



e-» = Ary^-4flg^cosh[Ai(»7 -»»)], z = z - (3-25) 

As before we can choose zq = without loss of generality. Although the metric, eq. (3.4), 
appears to be singular in this limit, this is only a consequence of an inconvenient choice for 
the time coordinate. If we instead convert to 'proper' time, r, defined by dr = t c / 2 dt, then 
r oc t l+c / 2 and the metric of eq. (|3.4p has a smooth large-c limit: 



ds 2 = _ e 2^) dr 2 + r 2 J e 2a(„) ff . . da ,i dx i + e Mv) d7] 2 + e »W de 2 . (3 . 26) 

We now convert back to 'conformal' time, t, using dr oc rdt, then r ~ e Hot ^ 2 , and 



2 _ c H t 



e 2w(v) dt 2 + e 2a(r,) § ^i^j + e 2v( V ) dr) 2 + ^6(^^21 ? (3 37) 



ds z = e 

As 7] -» ±00, e^"^ ~ (nq)' 1 ^^- - AH 2 e ±Xl ^-m) . Converting to proper radius dp = 
e v ^drj and identifying with the conical deficit form dp 2 + (1 — 5) 2 p 2 d# 2 we infer the deficit 
angles 

/V /^2 

<5 ± = 1 - Ai-^ °_ e ±Xi{m-V2) (3_2 8 ) 

Kg 

It is clear that by choosing any two of Ai, rji, rj2 and Hq appropriately, we can match this 
solution onto two conical branes of arbitrary tensions. 

Since in conformal time, all the 'scale factors' are exponential, this remains true in 4D 
Einstein frame. Such a solution arises when the effective 4D equation of state is w = —1/3, 
precisely at the transition from acceleration to deceleration. By itself this behaviour cannot 
be responsible for late time dark energy in the form of quintessence. 

Special case c = —2 

Option (ii) of ( 3.22j ) makes the choice c = —2, since in this case the (trj) Einstein equation 
implies the strong statement that £' = \{w — a) = everywhere throughout the bulk. For 
this special case the 6d geometry is everywhere maximally symmetric in the noncompact 4 
dimensions, taking the form 

ds 2 = e 2w ^ds 2 s4 + e^drj 2 + e 2b ^d0 2 , (3.29) 

where dS 4 is 4D de Sitter metric. These are just the de Sitter geometries considered in p3|| . 
Notice that there are no de Sitter solutions where both branes are conical. 
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3.3 Asymptotic forms 

We now re-examine the near-brane behaviour of these scaling solutions in order to connect 
their properties to those of the source branes. If we repeat the analysis of asymptotic forms 
given earlier for static solutions for the scaling solutions, with near-brane asymptotic form 
assumed to be given by eq. (|3.4|) , we find two changes relative to the static case. The simplest 
change is simply that the assumed time dependence implies that the coefficients c w , c a etc. 
now depend explicitly on t, with 

c 2 (t) oc c 2 a (t) oc t c , c 2 9 (t) oc t 2+c and cj(i) oc t" 2-c . (3.30) 

The second change is to do with the relationship amongst the powers a, (3, 7, uj and 
p which is dictated by the bulk equations. For instance, in the static case these equations 
required the powers to be related to one another by the Kasner-like conditions w + 3a + /5 = 
uj 2 + 2> a 2 + (3 2 + p 2 = 1. These conditions also apply for the scaling solutions, unchanged by 
c and Hq because the relevant terms in the field equations are subdominant in powers of p. 
However, in the scaling case there is also a new constraint, eq. ( |3.8| ), coming from the (trj) 
Einstein equation, which implies the following for the asymptotic powers: 

(2c + l)w + 3a + (2 + c)p = 0. (3.31) 

3.4 Generalised scaling solutions 

The previous scaling solutions admit a straightforward generalisation by allowing for the two 
compact dimensions to scale differently with time. The new ansatz is 

ds 2 = (H t) c \-e 2wM dt 2 + e 2a ^ 5ij dx l dx j ] + (H t) 2+C \e 2v ^d V 2 + (H t) 2s e 2b ^ d0 2 ] , (3.32) 

with 

, e v>(v) 

A e = {H t) s a e {r } ) and = — — ^— . (3.33) 

The technique for demonstrating that this is a solution follows as before. First we can show 
that all the powers of t drop out and the equations reduce to a system for the radial profiles 
alone. As before it is helpful to perform the redefinition ( 3.14 ) to the x,y,z,£ variables (we 
shall choose N = 1). In this case the (rjt) constraint implies 

f ' = — — (2e- 2x K 2 sa e a' e + 2sx' + (4 + 2c + s)z') . (3.34) 

b{2c + s) 

This can be solved for £ at least formally. In this case, we do not find a simple action for the 
system, but the remaining equations of motion can be expressed as 

- 2a' e x' - (1 + 2c)sH 2 e 2y+z '^a e = 
x „ _ e _ 6? _ 2x (^x+2y+z H 2 s{ y + 2c + s)+ e 2 y +z H 2 K 2 s 2 a 2 - e^K 2 {a' e ) 2 ) = 

y" - ±- e -<%-2*+2y (-2 e 2 *(2e 6 V - e z H 2 K 2 (2 + 2c 2 + s + s 2 + 2c(l + *))) + e z H 2 K 4 s 2 al) = 



z „ _ -^2 X+ 2y+z H 2 ( e 2. (6c _ g(g _ 3)) _ = ^ (3 35) 
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along with the second constraint 

3 (2c + s) 2 [(x') 2 + (z 1 ) 2 - {y') 2 + e- 2x K 2 {a' d ) 2 - ^e 2 ^ 

+e 2y+z-6£ H 2 ( 2(2 + 2( ,2 + ^ + ^ + ^ + 



+ 



(4 + 2c + s)z' + 2sx' + e~ 2x SK (a 2 ,)' 



0. 



(3.36) 



In this case it is necessary to resort to numerics to make further progress, nevertheless we 
have demonstrated that such a class of solutions exists. 

Scaling symmetry 

True to its name, the generalized scaling solution admits a scaling symmetry such that under 
t — > At, x — ► Ax and 9 — > \~ s 6, the metric transforms as g M]V — > X 2+c g MN . Infinitesimally 
this corresponds to the existence of a vector V 

d 



d d d d 
V = t 7^ + !~ h a^3^ — 

Ot OX\ OX2 OX3 



s9 



86 '' 



which satisfies the homethetic Killing vector condition [p7|1 

C-vQmn = —{D M V N + D N V M ) = —ag MN , 



(3.37) 



(3.38) 



for constant a = 2 + c. Spacetimes admitting such a vector are also known as self-similar. 
This vector is only globally well defined for s = since for nonzero s the interval 9 E [0, 2-7r] 
is not mapped onto itself. The dilaton and gauge field transform as eigenfunctions of this 
Killing vector 



Cve* 
C V A 



{2 + c)e\ 
(2 + c)A. 



(3.39) 



4. Noncompact conical solutions 

In this section we shall present some special time-dependent solutions that describe conical 
branes. Unlike the previous solutions, they describe a single brane in an uncompactified space. 
Nevertheless it is reasonable to expect that they could describe the near brane geometry of 
a more general time-dependent two-brane compact solution which is not described by the 
metrics already given. 

Our ansatz for the metric is motivated by the form of the scaling solutions: 

ds 2 = -a 2 (t)W 8 (t)N 2 (t)dt 2 + W 2 (t) (e^dx 2 + e~ 3 ^dr 2 ) + r 2 a 2 (t)d# 2 , (4.1) 

along with <j) = <f>{t) and gauge potential Ag = r 2 Q(t). Satisfying the (r,t) Einstein equation 
requires Q(t) = Q = constant and £ = (3y — 2z)/12 + £o- This is equivalent to saying 
W 2 e~ 3 t = d 2 e~ 3 ^ , and so we see that assuming we have normalized 9 to lie on the interval 
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[0, 27r], then the metric describes a conical brane with deficit angle 2n{\ — e 3 ^ ). It is useful 
to make the change of variables 

In a = (x - y + z)/4 (4.2) 

InW = (2x + y)/8 (4.3) 

(p = -{x + y + z)/2. (4.4) 

As usual the remaining equations of motion and constraint can be derived from the action 



S= dt 



N- 1 ( 6i 2 - 2>y - 2k 2 ) + N ( 24e 2 * 4 + 24 e 2£+2 ^V 



(4.5) 



This action looks qualitatively similar to the earlier actions, except for the signs of the kinetic 
terms. 

5. Brane and bulk dynamics 

In this section we relate the asymptotic form of the above solutions to brane properties, and 
use the result to argue that the solutions capture the late-time evolution of a pair of brane 
sources with nontrivial equations of state. 

5.1 Matching to brane properties 

It is possible to make a general statement of how brane properties dictate the asymptotic forms 
of the bulk fields in the near-brane limit, at least for branes for which gravity contributes 
negligibly to the total brane stress energy. This section makes this statement explicit for the 
6D case, following arguments presented in refs. S, S3]. 

A charged aside 

Before doing so, it is worth briefly pausing to develop some intuition from the analogous 
problem in electromagnetism. Consider for this purpose the electrostatic potential, y(r), 
generated by a collection of point charges, Qi, situated at various positions, rj(t), within 3 
spatial dimensions. In this case we know that some features of the resulting field are governed 
purely by the properties of individual charges, while others depend on the overall configuration 
of all of the charges. 

Typically the field very near the source charges depends purely on the properties of the 
nearby source, with the asymptotic behaviour having Coulomb form, <p(r) — ► Qi/\r — r$|, as 
|r— rj| — > 0. This form is ultimately dictated by Gauss' Law, which constrains the local electric 
field, E = — V</?, in terms of the local charge distribution. On the other hand, whether or not 
a given charge configuration is static — i.e. whether r-j = — is a type of question which 
cannot be purely determined using only near-charge properties, since it requires knowledge 
of the global positioning of all of the charges. 

This same kind of distinction arises also in the case of interest here: the gravitational 
fields sourced by a collection of branes. Again some features of the bulk geometry near the 
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brane are completely dictated by the physical properties of the brane, while others — most 
notably the time-dependence of the geometry — depend on the complete configuration of 
branes which are present. The next few sections identify which of the properties of the bulk 
fields are which. 



Thick branes and effective currents 

In this section we use the bulk field equations to show that local brane properties determine 
the near-brane form for radial derivatives (in the Gaussian-normal coordinates of § p.3| ) of 
the dilaton and gauge potential, d p <p and F„ a , as well as various combinations of the metric 
g a b and its extrinsic curvature, K a b = \ d p g a i,. Other quantities depend on the properties of 
all of the branes which source the solution, and so cannot be inferred purely from the local 
properties of nearby branes. 

To establish these points in the Appendix we review the arguments of ref. || (see also 
|32f ) and imagine regarding the brane source to be a 'thick' brane, which physically 
extends over a small proper distance, < p < I. Within this region we understand that the 
microphysical brane structure modifies the bulk equations, eqs. ( |1.2| ), to include new sources 
which are present only for p < £ and whose presence acts to smooth out the interior geometry 
at p = 0. Use of the field equations allows one to relate some of the near-brane properties of 
the external bulk fields to particular averages of these sources over the thick branes. 

If in particular these new sources depend only on the radial coordinate, involve only weak 
gravitational fields and the external bulk fields satisfy the asymptotic near-brane power-law 
behaviour of eqs. ( 2.11| ), then the constants in this asymptotic form are related to the brane's 
energy density, e = — t*t, pressures in the three noncompact and 6 directions, Pi = t l i and 
Pe = t dilaton 'charge', a, and Maxwell current, j 6 , by simple expressions: 



k 2 e 



K 2 p 
K 2 t e , 



2ir 
2tt 



2tt 

$' 



2tt 



and 



c e (uj + 2a + fiiHxCf- 1 - 1 
27Tc e (io + 3a)(Hi if' 1 = 2ttc (1 - ^{H^f' 1 

f = 27rc e c f (H 1 £y+P~ 2 . 



C0(1 -«)(#!<) 



13-1 



(5.1) 



These expressions, whose detailed derivation is given in the Appendix, identify which features 
of the near-brane bulk solutions are governed purely by the local properties of the brane, 
reproducing standard results in the case of a conical singularity with defect angle 2tt5 in the 
special case a = to = p = 0, (3 = 1 and eg = 1 — S. All of the remaining near-brane bulk- 
field properties — including in particular the time-dependence of the 4D metric — cannot be 
similarly determined purely from local information, and so depend in detail on the properties 
of both of the source branes. 

Notice that the above relations impose relations amongst the stress-energy components i a b 
once they are combined with the bulk field equations. In particular, the relations w+3 a+f3 = 
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1 and uj 2 + 3 a 2 + (3 2 + p 2 = 1 imply the components of the stress energy tensor must satisfy 



t a a = t l t + t' l i + t 6 e = — ;r Cg — 1 

i 

8vr 2 



(i^) 2 + £(^) 2 + (^) 2 + ^ = 2(<* " I) 2 + - /?) 



(5.2) 



where = cg(H£)^ 1 . For the simplest case of time-dependent scaling solutions discussed 
in 



3.1 , these must be supplemented by the additional relation following from eq. ( 3.31 ): 



(2 c + l)t\ + ^t i i -(2 + c)a = ^(C e - l)(c + 2) . 



(5.3) 



This equation is equivalent to the conservation of stress-energy on the brane. In general, in 
the presence of a dilaton coupling and current, the brane stress energy is not independently 
conserved but can couple with the bulk fields. 

Stress-energies for the scaling solutions 

We are now in a position to calculate the brane stress-energies that give rise to a specific 



scaling solution. For simplicity we shall consider the solutions of §3^. Since the near-brane 
limit corresponds to rj — > ±oo, it is sufficient to consider the asymptotic form of the scaling 
solutions, for which the metric components behave as exponentials. After a simple coordinate 
transformation, we can put the metric in the Kasner form 

ds 2 = (^) c [-/ t ± (^i^) 2w± dt 2 +/ a ± ( J H" 1 r) 2a± ^dxMx J 1+( J ff i) 2+c [dr 2 + f£(H ir ) 2P± - 2 r 2 d0 2 

(5.4)" 

Here the constants f t , and /« are determined by integration of the equations of motion. 
The Kasner exponents are given by formulae similar to the static solutions 



± 



a± 



ih 



Ai 



±5Af 
zfcA^ + 3Ai T 2Af 



±5A± 



Ai =f 2A^ 



and 



P± 



±Xf 



±5A± 



2(±A± 



Ai t 2A^ ' 
\iT2\f) 



±5Xf 



Ai T 2A 3 



(5.5) 
(5.6) 



where the A's are defined in §3.1. To directly compare with the boundary conditions, we must 
put the metric in GN coordinates. This can be achieved approximately with the transforma- 
tion r = (l?Qi) -1-c / 2 p, a t the cost of introducing new 0(p 2 ) contributions to the di 2 and didp 
components of the metric. These new contributions are negligible compared to those listed 
in eq. ( |5.4| ) in the p — > limit provided only that u < 1, which is in practice not a restrictive 
assumption. The result is a metric which takes the form 



ds 2 



-ft (Hot) 



c-w±(2+c) 



(H lP ) 2 ^dt 2 + ft{Hot) c - a±{2+c) (H lP ) 2a ±5 l3 dx l dxi 



+dp 2 + /±(F o i) (1 ^ ±)(2+c W) 2/3± d0 2 



(5.7) 
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and the dilaton and form-field are 



F P e = /±( J ff t)~ (2+c)(1+p±/2) (^iP) P±_1 



(5- 
(5- 



We note that it is important to transform to GN coordinates in order to make contact with 
the time-dependence of the brane geometry. In particular, whilst the coordinate position of 
the edge of the brane at p = £ is constant in GN coordinates, the f-dependence of the change 
of variables between r and p implies that it is instead located along a curve C(r,t) = I in 
another coordinate system. 

Finally we can read off the time-dependent brane stress-energies, dilaton coupling and 
current, 



k 2 e ± 



2vr [l - (1 - ^f^Hot^-^^HH^- 1 , 
-2n [l - (1 - a ± )ft{H t)^^){H^- 1 ] , 
n 2 pj = Ml - /3±)/, ± (^ot) (1 - /3±)(1+c/2) (^i^) /3± - 1 , 

3± = 27r/ e ± / / ± (i7 t) _(3+p± ^ ±)(1+c/2) (^^) /3±+p± • (5.10) 
In the special case of a conical singularity we have (3 = 1 and so a = to = p = 0, giving: 



K 2 e 



-n 2 pi 



2tt [1 - f e ] 



K 2 p e 



n 2 a 



0. 



2irf e ff(H t) 



-(2+c) 



(^1 



(5.11) 



The validity of the expressions ( 5.1C| ) rests on the weak- field assumption used in deriving 
the boundary conditions in the Appendix. Taking them literally for all time we would infer 
that the brane energy becomes negative at some point in the past or the future, something 
which is clearly unphysical. This arises because the effective deficit angle becomes negative. 
It seems likely that this is an artifact of our prescription for the boundary conditions, or 
that the bulk solutions are no longer valid in this regime. To resolve this it will be necessary 
to build smooth, defect-like models of the branes and analyze the relationship between the 
boundary conditions and the bulk solutions. We feel that nevertheless, the scaling solutions 
we have found here will at least represent a consistent description of the bulk in one or more 
asymptotic regimes. 



Delta-function sources and conical singularities 

It is a common practice to represent the brane sources in terms of a point-like delta-function 
in the two extra dimensions, with an action of the form 



S b 



d 6 xV=7 !(<]>) 5 2 (x), 



(5.12) 
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where 7^ = g MN d^x M d u x N denotes the induced metric on the brane. In this case direct 
variation of the brane action with respect to bulk fields would lead to 

T (b) ~ *(6) <V 6 " ' J (&) " Hb) 5 » and P {b) - , I 5 - 13 ) 

with 

*SS = -hi£r = ~ m ^ ' a w = — ^ = //((/,) ' (5 - 14) 

and f {b) = -5S b /5A u = 0. 

However, it is important to realize that such an assumption requires that the bulk fields 
remain regular at the brane position, p —* 0, so that it makes sense to evaluate the bulk fields 
there. The above expressions show that this is not generic for branes having co-dimension 
> 2, since it requires the powers a, /3, u and p all to be non-negative. If the brane couples 
to the Maxwell field then the condition on p strengthens to p > 1, which is consistent with 
the bulk field equation a 2 + (3 2 + u> 2 + p 2 = 1 only if p = 1 and a = j3 = to = 0. Otherwise 
the representation of the brane in terms of a (^-function source can be overly restrictive. 
More generally one can instead describe the low-energy bulk dynamics by excising the brane 
positions from the bulk geometry, and describing the brane sources in terms of boundary 
dynamics on the resulting co-dimension one boundary of the excised region |^9[ . 

5.2 Parameter counting 

^provides an 8-parameter family of static solutions, with independent parameters A2, A3, A4, 
Q, VI1V2, zq and £o- ^1 can be fixed by means of the constraint equation A| = \f +A|+12 A|. Of 
these, two (say r/2 and £0) can be removed by appropriately choosing units in the noncompact 
four dimensions, since this can be used to set a and w to any convenient value at any one 
position. (Notice that 6(770) cannot be similarly set to vanish in this way without changing 
the periodicity condition 6 ~ 6 + 2tt. We also do not use the classical scale invariance of 
the bulk field equations to remove zq by allowing <j) to be set to any desired value at any one 
point, i] = rjQ. We do not do so because this scale invariance may be broken by the brane-bulk 
couplings, and so may not be consistent with the boundary conditions given in the previous 
section.) This leaves a total of 8 — 2 = 6 physical parameters in the static solutions presented 
above. 

A similar counting applies to the generalised scaling solutions, whose time dependence 
introduces three more parameters: c, s and Hq.® Integrating the remaining equations to 
obtain the field profiles as functions of rj may be expected to introduce only 5 more constants 
(rather than 6, due to the additional Bianchi-identity related constraint relative to the static 
case) into the general scaling solution as well, leading us to expect there to be a total of 
5 + 3 = 8 parameters in these solutions. 

We may now count the number of parameters which should be expected of a general 
bulk configuration which is sourced by two branes. For the symmetries of interest each brane 



"We note that Hi is not an independent parameter. 
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is characterized by the 5 quantities e, pi, pg, j e and a. With two source branes there are 
10 such parameters in total. To this should also be added the integer which measures the 
total magnetic flux through the extra dimensions, so long as we also include the topological 
constraint [10| which relates this integer to the magnetic currents on each brane. This leaves 
us with a total of 10 independent parameters describing the physics which sources the bulk 
fields, one of which is quantized to be an integer. This represents 2 more parameters than our 
general time-dependent solutions have available to accomodate. We believe this is because 
the solutions we consider are not singular enough to contain direct couplings between the 
branes and the Maxwell fields, and so j 8 = on both branes. 

Scaling solutions are often attractor solutions towards which general time-dependent 
configurations tend after the passage of any initial transients. If this is also so for the 6D 
supergravity field equations, we would be led to the following attractive picture. It has long 
been known | 10| , 22, 12] that static (and maximally-symmetric, but curved ]23|| ) solutions can 
only exist if the properties of the two source branes are appropriately adjusted relative to one 
another. But it has been unknown what happens to the bulk geometry in the generic case 
where such adjustments are not made, although it has been suspected that these would pro- 
duce time-dependent bulk configurations. Based on the above considerations, in the generic 
case we expect that the bulk is indeed time-dependent, and in particular this time-dependence 
approaches one of the scaling solutions given here (once transients pass) at late times. 



6. Conclusions 



Six dimensional supergravity provides a fascinating laboratory for investigating the issues 
which underly the cosmological constant problem, largely because six dimensions is both 
simple enough to allow the construction of explicit solutions, yet rich enough to exhibit an 
interesting variety of properties. In particular, it provides the simplest setting within which 
a collection of positive-tension branes can combine to produce vanishing 4D curvature. This 
makes it a very fruitful arena in which to explore how natural are the choices which must be 
made in order to ensure acceptably flat 4D worlds. 

Our main result in this paper is to provide a new class of static and time-dependent solu- 
tions to the full field equations of gauged chiral 6D supergravity. For both classes of solutions 
there are 4 warped, noncompact dimensions (which need not be maximally symmetric), and 
2 curved compact dimensions which are taken to be axially symmetric. The time-dependent 
solutions are of the self-similar scaling type, and as such are plausibly attractors which ro- 
bustly describe the late-time limit of the time-dependent evolution generated by fairly generic 
initial conditions. 

The solutions we find diverge at (at most) two positions within the extra dimensions, 
which we interpret to be the positions of two space-filling 3-branes. This interpretation is 
supported by examining the asymptotic behaviour near the branes, which has a relatively 
simple interpretation in terms of brane stress-energy and dilaton coupling, at least in the 
limit that the gravitational fields involved in the brane structure are very weak. 
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The picture which emerges of brane-bulk dynamics is this. Much like for a collection 
of static electric charges within a compact space, the fields produced by a generic set of 
sources is time dependent. It is nonetheless possible to achieve static solutions provided the 
properties of the two branes involved are appropriately adjusted relative to one another. In the 
supersymmetric case the resulting static solutions are also marginally stable (classically) [30| . 
[HJ], with the marginal perturbations being along a flat direction whose presence is guaranteed 
by a classical scale invariance of the 6D supergravity equations. We believe our scaling 
solutions to be the natural endpoint for motions which start out along these flat directions. 

What do these results imply for the 6D self-tuning mechanism for the cosmological con- 
stant? It must be emphasized that they do not at all address the central issue of technical 
naturalness — i.e. the question of whether the choices required for flat (or slowly- varying) 
cosmologies are stable against renormalization. What they do address is the ancillary issue 
of initial conditions, which arises in any theory for which the Dark Energy density is evolving 
in time. They do so by identifying how special the solutions are which can give acceptable 
4D cosmologies. 



Although such acceptable solutions do exist 26 ] , their existence requires both an accept- 
ably shallow potential and logarithmic corrections to this potential (such as can be generated 
by Casimir energies in the 2 extra dimensions [14|). However the slow roll due to potentials 
like these coming from a quantum origin can only be relevant if they are not dominated by the 
fast motion driven by larger classical forces, and the motion described by the generic solutions 
obtained here would be much too fast to provide an acceptable 4D cosmology. So any descrip- 
tion of the Dark Energy in terms of 6D dynamics along these lines must presuppose initial 
conditions which exclude the generic motion found here, and so rely on the universe being 
prepared with branes whose properties are sufficiently well adjusted relative to one another 
to give close to static classical dynamics. One might hope that such special initial conditions 
might have their explanation in the same way as did the initial condition problems of the 
standard Hot Big Bang cosmology: that is, in terms of the (possibly inflationary) dynamics 
of still-earlier epochs of the universe. 



Acknowledgements 

AJT is supported in part by US Department of Energy Grant DE-FG02-91ER40671. CB 
would like to thank the Banff International Research Station for its kind hospitality while 
this paper was being finished. CB and DH acknowledge support from the Natural Sciences 
and Engineering Research Council of Canada, with additional support for CB coming from 
the Killam Foundation and McMaster University. DH acknowledges support from les fonds 
Quebecois de la recherche sur la nature et les technologies. CdR is funded by a grant from the 
Swiss National Science Foundation. Research at Perimeter Institute is supported in part by 
the Government of Canada through NSERC and by the Province of Ontario through MEDT. 



- 21 - 



A. Brane bulk matching 



In this Appendix we summarize the connection between the near-brane asymptotic behaviour 
of the bulk solutions and some of the physical properties of these branes, following the argu- 
ments made in refs. [22, §[|. To this end we choose Gaussian-normal coordinates as described 
in C 



Our goal is to use the bulk field equations to show that local brane properties determine 



the near-brane form of the radial derivatives of the bulk fields: d„ 



F pa and K ab 



\ d P g a b- 

Other quantities, such as the time dependence of the near-brane metric, depend on the 
properties of global properties involving all of the branes which source the solution, and 
so cannot be inferred purely from the properties of the brane at p = 0. 

To establish these points we follow refs. 12, @, 32] and imagine regarding the brane 



source to be a 'thick' brane, which physically extends over a small proper distance, < p < i. 
Within this region we understand that the microphysical brane structure modifies the bulk 



equations, eqs. (1.2), to include new sources in such a way as to ensure that all quantities 
remain smooth as p — > 0. We quantify this smoothness by assuming that the resulting bulk 
fields obey the conditions = 0(p 2 ), Kgg = 2p, F ap = O(p) and d p <f> = O(p) as p — > 0. 
We define the new sources which are required in this way from the equations 



Gain ~t~ K T*. 



D M F 



□ « 



K 2 P 



o, 



(A.l) 

where G MN is the Einstein tensor. The source currents J N = J { N B) + J® and P = P^ + i>W 
include both the contributions of the bulk fields, l (P>)\ as inferred from eqs. (1.2) plus new 
brane contributions, '(6)', which vanish for p > I. The stress energy receives similar kinds of 
contributions, T MN = + T$ N . 

Now comes the main argument. Consider for simplicity a purely radial source profile 
within the thick brane, with 

ds 2 = g^(p) dx^dx" + g mn (p) dx m dx n = g^(p) dxW + e 2b ^ d9 2 + dp 2 , (A.2) 

and so on for the other bulk fields. The point of the definitions given above is that the field 
equations allow the sources to be written as total derivatives, allowing their averages over the 
thick-brane volume 7 to give expressions depending only on the boundary values at p = and 
p = i. This is simplest to see for the dilaton, for which 



K 2 a 



p<1 

2-k 



d9 



^K 2 P = 


f d 2 x v 






Jp<£ 




Jgi s M d M (j) 


p=i 

= 2vr 






p=0 





p=l 



(A.3) 



where s M is the outward-pointing unit normal to surfaces of fixed p and the last equality uses 



the above-mentioned requirement that d p <j) must vanish at p = 0. Eq. (A.3) gives an explicit 



relation between d n 



and purely local brane properties. 



7 Multipole generalizations can be similarly denned by appropriately weighting the integrands by functions 



of 0. 
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A similar argument may be made for the current J M . Averaging over the thick brane 
volume gives 



f = I d 2 x^J a =[ d 2 x^D M F 

J 0<i J 0<l 

fg~2s M F Ma 



2tt 



d9 



p<i 
p=i 



2ir 



p=Q 



b ppa 



(A.4) 



with the contribution from p = again vanishing due to the condition that F pN must vanish 
there. We see from this how the near-brane limit of F pa relates to local brane properties. 

A similar result for the energy density is more difficult to define in general, however 
progress is possible for purely radial profiles if the gravitational binding energy of the thick 
brane is negligible. In this case, averaged stress energy components can be defined by 



n 2 t a b = [ d 2 x^K 2 T a b = - [ d 2 x^G a b 
J p<i Jp<e 



(A.5) 



ip<e J p<t 

It is useful at this point to specialize to the particular metric, ds 2 = —e 2w ^ dt 2 +e 2a ^ 5ij dxMa>' + 
e 2b (p) (\Q 2 _|_ dp 2 , for which ^fg^ = e b and the nonvanishing components of the Einstein tensor 
are 



-G\ = 3 a" + b" + 6 (a') 2 + 3 a'b' + (b'f 

-G i j = \w" + 2 a" + b" + 2 w'a! + 3 (a') 2 + 2 a'b' + {w') 2 + w'b' + {b'f 
-G% = w" + 3 a" + 3 w'a' + (w') 2 + 6 (a') 2 , 



(A.6) 



where ' = d/dp. The main point now is this: although the above integrals depend on the 
detailed metric profiles, the result simplifies in the limit of weak gravitational fields for which 
the expression may be linearized in the derivatives w', a' and b , where e 2b = p 2 e 2b . In this 
case the integrals may be performed explicitly, leading to the expressions 



2 t 

K t t ~ 27T 



K t j ~ 27T 



2 6 

k t $ ~ 2tt 



e b (3a' + b') 



p=0 



e b (w' + 2 a' + b') 



p=o' 



e b (w' + 3a') 



p=0 



(A.7) 



We now specialize the above expressions by evaluating them using the smooth limit at 
p = and the asymptotic form of the bulk metrics near p = I: 

ds 2 ~ -[c w {HpY\ 2 dt 2 + [c a {Hp) a ] 2 5ij dx { dx j + dp 2 + [c e (Hpf- 1 ] 2 p 2 de 2 
F pe ~c f {H P y and ~ c^Hpf . (A.8) 

Here H is an arbitrary dimensionful scale, while u, a, (3, p, c w , c a , c$, Cf and are constants, 
for which the bulk field equations imply the conditions uj + 3a + (3 = u 2 + 3a 2 + (5 2 +p 2 = 1 and 
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7 = V ~ 1- For this metric the relevant derivatives are w'(£) = uj/£, a!(l) = a/£, b'(£) = (5/1 
and cf>'{t) = p/£, while the 2D volume element is ^fg^ = eg (Hp) 13 . 

The relevant quantities at p = are e b w'\o = e b a'\o = and e b b'\o = 1, and so up to 
contributions that are sub-leading for small £ we find 



2tt 
2vr 



ceiSa + p^mf- 1 -1 
c 9 (w + 2a + /?)(«) /3 - 1 



2tt 



(l-w)^)^ 1 -! 



2tt 



c e (l-a)(^) 



/3-1 



K 2 t% ps 2TTc e {uJ + ^a)(H£f- 1 
K 2 a = 27rc e p(H£f- 1 
and / = 2ircgCf{H£)~< + P- 1 . 



2tt c»(1 - (3){H£) 



(A.9) 



These are the expressions which are used in the main body of the text. Notice that in 
the case of a purely conical singularity having defect angle < 5 < 1 in the bulk geometry we 
have a = u = p = and (3 = 1 and eg = 1 — 5, in which case the above expressions simplify 
to the following results for the energy density, e, 3D pressure, pi, and off-brane pressure, pg: 



k 2 £ 



-27T(5 and n A pg = K 2 t 6 g . 



(A.10) 
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